Chebyshev tau-QZ algorithm methods for calculating spectra of hydrodynamic stability problems

The Chebyshev tau method is examined in detail for a variety of eigenvalue problems arising in hydrodynamic stability studies, particularly those of Orr-Sommerfeld type. We concentrate on determining the whole of the top end of the spectrum in parameter ranges beyond those often explored. The method employing a Chebyshev representation of the fourth derivative operator, D 4, is compared with those involving the second and first derivative operators, D 2 and D, respectively. The latter two representations require use of the QZ algorithm in the resolution of the singular generalised matrix eigenvalue problem which arises. Physical problems explored are those of Poiseuille flow, Couette flow, pressure gradient driven circular pipe flow, and Couette and Poiseuille problems for two viscous, immiscible fluids, one overlying the other.